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£N1 ' Abstract 

> ' 

We present a new method to measure the transverse spin correlation 
in the H — » t + t~ decay. The method has been devised to be insensitive 
f"^ ■ to the beamstrahlung which affects the definition of the beam energy 

I/-) ' at a linear collider. In the case of two r —* ir v t {vt) decays, using 

f*"^ . the anticipated detector performance of the TESLA project, we get a 

' promising estimation of the error expected on the measurement of a CP 

violating phase. 



1 Introduction 



The possibility to determine the CP properties of a light Higgs boson through 
the spin correlations in its H — > t + t~ decay has been often considered [1-9]. 
The principle is simple. Let ± denote the projection of the spins of the r's in 
. their respective rest frames on a z-axis oriented in the direction of the r _ line 

5— i ' of flight for the r + and opposite to the r + line of flight for the r~ . The t + t~ 

spin state 

^[| + -> + e«|-+>] (1) 

is a CP= +1 state for £ = 0, a CP= —1 state for £ = it, and a mixed CP state 
otherwise. Such a state is produced by the decay of a CP=+1 Higgs with a 
coupling [5, 9] 

gf(cosif) + i sin ipj^rH . (2) 

In this case, neglecting 0{m^/m 2 H ), the phase is £ = 2ip. 
The spin correlations for the state (1) are: 

C zz = -1, C xx = C V y = cos£ = cos2ip, C xy = -C yx = sin £ = sin2^ . (3) 
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The way to measure ip (£) is transparent in the case of two — > tt^P t (u t ) 
decays. Using (1) and the r decay amplitudes, one gets the correlated decay 
distribution 

W 3 (cos # + , cos #~ , Acp) = -!- [l + cos 6>+ cos 6»~ - sin 6 + sin 6>~ cos( Atp - 2ip)] , 

(4) 

where 9^ is the polar angle in the rest frame between the tt^ direction (ir) 
and the above defined z-axis. The relative azimuthal angle Aip = (p + — (p~ is 
the angle between the two planes defined by the t~ direction and the ir + 
direction respectively in the Higgs rest frame. The distribution of the azimuthal 
angle is obtained by integrating out the polar angles: 



Wl (A v ) = j- 



ir 2 

1 - — cos(A^ - 2V) 



(5) 



The merits of the two distributions for the measurement of ip can be quantified 
by their sensitivities = l/cr^V^V, where is the error on ip expected from 
a maximum likelihood fit of the distribution for a sample of N events. The 
sensitivities measure the information per event on tp contained in the distribu- 
tions [10]; their computation is straightforward when an analytical or numerical 
expression of the distributions is known. They are 5^ = 1.15 for the distribu- 
tion (4) and — 0.92 for (5). The superiority of (4) is not very large and it 
decreases when experimental effects are introduced. For that reason and for the 
sake of simplicity, we will only consider the distribution (5) in the following. 

The best place to study a light Higgs at a linear e + e~ collider is the Hig- 
gsstrahlung process e + e~ — ► ZH. We assume that the mass of the Higgs is well 
measured by the analysis of its dominant decay modes. The four- momentum of 
the Higgs is determined by the measurement of the Z and therefore the four- 
momenta of the t + and the r~ can, in principle, be reconstructed in the case 
of two hadronic decays (Sec. 2). It is known [10] that under such circumstances 
the three main hadronic decay modes, — > 7r ± i? T (^ T ), — ► p il v T (y T ), and 
r ± — > afv T (v T ) have the same analysis power for the measurement of spin ef- 
fects. This is due to the fact [11,13] that, in the r rest frame, all the information 
on the spin is embodied in the distribution of a unit vector a, the polarization 
analyser, which is equal to 7r (pion direction) in the case of — ► 7r ± f' T (^ T ). 
This vector can be computed [12, 13] from the measured momenta in the case 
of the above decay modes. For example, if r ± — > p ± D T (i> T ) and p ± — > 7r ± 7r°, 
then 

a i =Af(2( q -p v )q i -(q-q)pi) , (6) 

where J\f is a normalization factor, p u = p T ± —p p ± is the four- momentum of the 
neutrino, and q = p^± — p^a is the difference of the four-momenta of the two 
pions. The distribution (4), where tt is replaced by a in the definition of the 
angles, contains all the available information on the spin correlation. The three 
decay modes (i.e., (55 %) 2 of the t + t~ pairs) can therefore be used in the same 
way for the measurement of the CP violating parameter tp [5] . 
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Unfortunately, it has been shown by Monte Carlo studies [6] that the stan- 
dard reconstruction (Sec. 2) of the r four-momenta is critically impaired by the 
effects of beamstrahlung. A new estimator has been proposed [7-9] , which is 
less sensitive to the quality of the r reconstruction but requires that both r + 
and t~ decay by the process — > p^D T (v T ). However, the ideal sensitivity 
of the new method, when all the four-momenta are exactly known, is 0.48 to 
be compared with 0.92 for the standard method. Besides it takes advantage of 
(25%) 2 only of the r + r~ pairs. There is therefore a possible improvement of 
the error on ^ by a factor of at least four, if a method of reconstruction of the 
r's insensitive to beamstrahlung is found. 

The aim of the present paper is to devise such a method. To check the 
sensitivity to beamstrahlung, we use a simple Monte Carlo procedure. First 
the energies of the e ± after beamstrahlung are generated using the program 
circe [14] with the parameters of the TESLA project and the cross section for 
Higgsstrahlung in the case of a CP=+1 Higgs [4]; next the production and decay 
angles of the Z are generated according to their correlated distribution under 
the same hypothesis [4]. The Higgs decay angles are generated according to 
an isotropic distribution. For the correlated decay of the r's we consider three 
cases: two decays into ttv (tttt), a decay into tti> and a decay into pv (/wr), and 
two decays into pv (pp). For each case, the azimuthal angles and the cosines 
of the polar angles in the decays of r's and p's are generated uniformly. From 
Eq. 1 and the decay amplitudes for r* — > 7r ± i> r (^ r ) and r^ — > p ± v T {v T ) — > 
7r ± 7r°^ r (v T ), the correlated decay probabilities are computed for diverse values 
of tp and the events are accepted or rejected accordingly. The four-momenta of 
all the particles are then computed and, finally, the decay lengths of the r's are 
generated. 

The Higgs mass is assumed to be 120 GeV and three values of the total energy 
are considered: </s = 230 GeV, V* = 350 GeV and y/s = 500 GeV. The natural 
energy for a detailed study of the Higgs is the energy of the largest cross section 
(i.e. near 230 GeV). The consideration of higher energies allows estimating 
the robustness of the method and possibly the effect of an underestimation of 
radiative effects. 



2 Reconstruction of the r's in the Higgs rest 
frame 

The reconstruction of the r's in the r + r~ rest frame for hadronic decays t ± — > 
h ± v T (v T ) has been known for a long time [15] and was used at LEP to improve 
the measurement of the r polarization [16]. Its principle is sketched in Fig. 1. 
Both r's have the same energy = and the energies and momenta of 

the hadrons are measured. The angle 8^ between the direction of a r ± displayed 
by the unit vector f ± and the direction of the hadron fi^ is therefore fixed: 
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Figure 1: Reconstruction of the r direction in the Higgs rest frame. 



The direction of the r~ , f = t~ is on the intersection between the cone around 
h~ with angle 5~ and the cone around — h + with angle S + . 

In general there are two solutions f 1 / 2 = n y ± n±, where the vectors n y and 
nj_, shown in Fig. 1, are given by 



cos S + — cos S cos 5 f , cos S — cos <5 cos S 



n« = -- 



h+ + r - y - /i- , (8) 



sin <5 



sin 2 (5 sin 2 5 

■ h + A h~ , with cos S = —h + ■ h~ 



The ambiguity can be resolved using the information from a vertex detector 
but, because the main vertex is known from the Z decay the use of the detector 
is rather different and simpler than at LEP [16, 17]. One needs only to make a 
chi-square test on the distances in the laboratory between the reconstructed r 
lines of flight and the trajectories of the charged pion's. 

The situation is degraded in the presence of beamstrahlung and/or other ex- 
perimental effects because the intersection of the two cones is no more granted. 
As a result the acceptance is reduced and the distribution of Atp strongly de- 
formed. The reconstruction is especially awkward in the case of two — ► 
n ± D T (i' T ) decays, because the polarization analyser is then a ± = 7r ± = h ± and 
thus Aip is the angle between the vectors A h~ and t~ A h + . From Eq. 8, 
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one gets 



1 - cos 2 Aip = (1 



sin S 



■ TV 



(9) 



sin o + sin d 

which shows that the two cones are tangent when | cos Aip\ — 1, feature that can 
also be observed in Fig. 1. As a consequence, even for yfs = 230 GeV, where the 
closeness to the threshold reduces the effect of beamstrahlung, the acceptance 
becomes very small when cos Aip is near ±1 and the distribution (5) can hardly 
be used to test the CP properties of the Higgs. 



3 Reconstruction of the r's in the laboratory 

Owing to the observation of the Z decay products, which allows the reconstruc- 
tion of the main vertex, one may envisage to perform the reconstruction of the 
r's in the laboratory frame. 

Let us assume that the t ± energies in the laboratory E^± are known. The 
angles between the directions of the r's and the hadrons are given by the 
relation (7), which reads here 



cos a 



(10) 



The reconstruction of the r ± direction is then very simple in the case of the 
— > -k ± v t {v t ) decay mode (Fig. 2-a). 





Figure 2: Reconstruction of the t direction in the laboratory frame: (a) for a 
r ± — > 7r ± i? T (^ T ) decay mode, (b) for a t ± — > 7r ± 7r°f' r (^ T ) decay mode. 



Denoting by I the unit vector of the perpendicular from the vertex to the 
pion trajectory and by 7r, the unit vector of the pion momentum, the vector f 
is given by 

f = cosa7r + sin a i . (11) 
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In the case of a — > 7r ± 7r°f' T (^ T ) decay mode, f is on the intersection of the 
plane defined by the vertex and the charged tt trajectory with the cone around 
h with angle a (Fig. 2-b). There are in general two solutions. The problem 
of the reconstruction of the two r's is therefore a problem with two unknowns: 
E^ + and E^_ , but even in the presence of beamstrahlung we still have three 
constraints: the conservation of the components of the momentum orthogonal 
to the beams and the equality of the t + t~ effective mass with the Higgs mass. 
They are sufficient to determine the two energies and resolve the ambiguities if 
needed. 

4 A simplified algorithm 

Implementing the last method by a fit would require a good knowledge of the 
errors and their correlations. This is not possible with our simple simulation. 
For that reason, we will use a new procedure, which combines elements of the 
approaches followed in the two previous sections and gives good results, without 
the intricacies of a fit. 

Taking pn = —pz, we start the reconstruction in the Higgs rest frame but 
use the approximation f = n||/|nj. This is always possible if cosc5 ± is replaced 
by 1 when it is found greater than 1. 

As the r energies in the Higgs rest frame are known: E T ± = mH/2, we can 
compute E^ ± and perform the r reconstruction in the laboratory. For that, 
we replace cos by 1 when it is found greater than 1 and, in the case of a 
r ± — > p ± v T (v T ) decay, use the projection of h on the (V,7r,t) plane to define 
f when the cone and the plane do not intersect. In the tttt channel, an event 
is rejected if both cosa + and cos a - are greater than 1. The ambiguities are 
resolved by choosing the solution with the smallest missing \p±\. We can now 
redefine the Higgs frame as the rest frame of the t+t~ pair and compute Aip. 

The reconstructed distributions for the three channels at an energy of ^/s = 
350 GeV are shown in Fig. 3 and their sensitivities to at the three considered 
energies are given in Table 1. 

Both the curves in Fig. 3 and the numbers in Table 1 include the effect of 
the small loss of acceptance due to the rejection of events in the tttt channel. 



Table 1: The sensitivities to ip of the reconstructed distributions when only 
beamstrahlung is taken into account. The effect of the small loss of acceptance 
in the tttt channel is included. 



Vi (GeV) 


Sensitivity (S^) 




TTTT TT P pp 


230 


0.92 0.88 0.83 


350 


0.91 0.73 0.66 


500 


0.88 0.64 0.55 
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Figure 3: The distributions of A<p: (a) at the generation level, (b), (c), and (d), 
reconstructed by the method of Sec. 4 for the three channels mr, pir, and pp at 
an energy of 350 GeV. Beamstrahlung effects only are taken into account. The 
histograms are normalized to the number of generated events and multiplied by 
the number of bins. With this normalization, the distribution (a) is 2ttWi(Aip). 
The full lines correspond to tp = 0, the dotted lines to i\) = n/8. 



For the irir channel, the sensitivity is nearly the ideal one (0.92) up to 
500 GeV. For the pit and pp channels, the sensitivities are slightly reduced at 
230 GeV and decrease more rapidly with the energy than for the irir channel. 
Two effects contribute to that. The first is the closing of the r decay angle 
when the hadron mass increases. The second is the imperfect resolution of the 
ambiguities. The second point can be improved because the criterion used for 
the choice of the solution is not optimal. Other constraints can be used, like 
the t + t~ effective mass and the positivity of the decay length. It should also 
be noted that the conservation of p± is an important but not vital point in the 
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method. For example, smearing the p± of the Z with a(p x ) = cr(p y ) = 1 GeV 
yields sensitivities S p7T = 0.66 and S pp = 0.52 at y/s = 350GeV. 

Finally, a few remarks are in order about the — > afv T {v T ) decay channel, 
which is not included in our unsophisticated simulation, for the reason that its 
description is more complex than a simple angular distribution. The af — ► 
tt^~ it® it® decay mode is reconstructed by the same method that the p ± , but 
the sensitivity will probably be worsened by the closing of the r decay angle. 
For the af — > 7r ± 7r + 7r _ decay mode, the determination of f in the laboratory 
is in principle possible from vertexing information only, consequently a good 
sensitivity can certainly be achieved by using this determination of f or by an 
adaptation of the method used for the p ± . 



5 A semi-realistic simulation of the mz channel 

To get a realistic estimation of the sensitivities, it is necessary to take into 
account the performance of the detector. We use for that the parameters of the 
TESLA project [18, 19]. Since the precision of the ir° measurement depends not 
only on the accuracy of the detector but also on the quality of the reconstruction 
algorithm, we consider here the tttt channel only. 

For the charged tracks an independent Gaussian smearing is performed on 
the five parameters: 6, <fi, l/p± and the two components of the impact param- 
eter. We use for the widths of the Gaussians the following values [19]: 

a(6) = a(<p) = 0.1 mrad, cr(l/p ± ) = 5 x 10" 5 GeV" 1 
and a(r(f>) = a(rz) = (4.2 ® 4.0/(psin 3/2 6>)) A™ • 

The energy of the jets is smeared accor ding to a(E)/E = 0.3/V^(GeV) . The 
position of the vertex is determined by the shape of the beam [18] for the x and 
y coordinates and by the charged decay products of the Z for the z coordinate. 
The smearing of these coordinates is done accordingly. 

The reconstructed distributions at an energy of 350 GeV, with all the exper- 
imental effects included in the simulation, are shown in Fig. 4 for the decays of 



Table 2: The sensitivities to tp of the reconstructed distributions for the itit 
channel, when all the experimental effects are taken into account. 



Vi (GeV) 


Sensitivity (5^) 




Z — > fi + Z — > qq 


230 


0.69 0.71 


350 


0.60 0.61 


500 


0.58 0.58 
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Figure 4: The reconstructed distributions of Aip for the irir channel at an energy 
of 350 GeV, when all the experimental effects are taken into account: (a) for 
the decay of the Z into (b) for its decays into two jets. The convention 

of normalization and the values of ip are the same as in Fig. 3. 



the Z both into and into qq. Their sensitivities to tp at the three consid- 

ered energies are given in Table 2. It is clear from these values that a large part 
of the sensitivity is retained. 

Because of the key role of the vertex detector in the reconstruction, we have 
done again the simulations with the vertex detector errors multiplied by two. 
The reduction of the sensitivities that results from the increased uncertainties 
is always smaller than 0.1. 



6 Conclusion 

We have studied the production of a light Higgs boson by the process of Hig- 
gsstrahlung and its subsequent decay into t + t~ , under the conditions of a linear 
collider. 

We have described a method, which by the joint use of kinematics and 
vertexing allows the measurement of the transverse spin correlations of the two 
r's. This method is not impaired by beamstrahlung and can be applied for the 
main hadronic decay modes of the r and most of the visible decay modes of the 
Z. 

In the case of two t ± — ► n ± D T (v T ) decays, a complete simulation of the detec- 
tor effects with the parameters of the TESLA project [19] has been performed. 
A realistic study of the reconstruction of the ir° is still to be done, neverthe- 
less it appears that a reasonable goal for the measurement of the phase tp that 
parametrizes a possible CP violation should be to use all the above mentioned 
final states and get a sensitivity better than 0.5, i.e., cr^ < 0.6 n/y/N evt .. 
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